Slowly rotating pulsars 
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In the present work we investigate one possible variation on the usual static pulsars: the inclusion 
of rotation. We use a formalism proposed by Hartle and Thorne to calculate the properties of 
rotating pulsars. All calculations were performed for zero temperature and fixed entropy equations 
of state. 
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I. INTRODUCTION 

Pulsars are believed to be the remnants of supernova 
explosions. They have masses 1-2M©, radii ~ 10 km, 
and a temperature of the order of 10 11 K at birth, cool- 
ing within a few days to about 10 10 K by emitting neutri- 
nos. Pulsars are normally known as neutron stars. Qual- 
itatively, a neutron star is analogous to a white dwarf 
star, with the pressure due to degenerate neutrons rather 
than degenerate electrons. The assumption that the neu- 
trons in a neutron star can be treated as an ideal gas is 
not well justified: the effect of the strong force needs 
to be taken into account by replacing the equation of 
state (EoS) for an ideal gas by a more realistic EoS. The 
composition of pulsars remains a source of speculation, 
with some of the possibilities being the presence of hy- 
perons [lL |sl H[ , a mixed phase of hyperons and quarks 
[1, IE Bui la]) a phase of deconfined quarks or pion and 
kaon condensates @. Another possibility would be that 
pulsars are, in fact, quark stars [111 ]. In conventional 
models, hadrons are assumed to be the true ground state 
of the strong interaction. However, it has been argued 

mis mil! 

Htjj that strange matter composed of de- 
confined u, d and s quarks is the true ground state of all 
matter. In the present work we have opted to use 
the term quark stars instead of strange stars be- 
cause in more refined nuclear structure models, 
as the Nambu-Jona-Lasinio model for instance, 
the s quark appears in a much smaller quantity 
than in the very naive MIT model, where it is 
con sidered in equal footing with u and d quarks 
| 111 ] . In the stellar modeling, the structure of the star 
depends on the assumed EoS, which is different in each of 
the above mentioned cases. An important distinction be- 
tween quark stars and conventional neutron stars is that 
the quark stars are self-bound by the strong interaction, 
whereas neutron stars are bound by gravity. This allows 
a quark star to rotate faster than would be possible for 
a neutron star [ll|, [l6| . 

Once an adequate EoS is chosen, it is used as in- 
put to the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tions fl7j | . which are derived from Einstein's equations 
in the Schwarzschild metric for a static, spherical star. 
Some of the stellar properties, as the radius, gravitational 



and baryonic masses, central energy densities, etc are ob- 
tained. These results are then tested against some of the 
constraints provided by astronomers and astrophysicists 
[l8l HU and some of the EoS are shown to be inappropri- 
ate for describing pulsars [!, IS ■ 

On the other hand, it is well known from the 
Doppler broadening of the pulsars spectral lines 
that they rotate. Some pulsars are observed to ro- 
tate with periods as small as 1.56 ms [10]. Newly 
born hot pulsars can rotate fast and hence un- 
dergo instabilities. The effect of a rotating star on 
spacetime is commonly known as the Lense-Thirring ef- 
fect, which refers to the dragging of local incrtial frames. 
Rotation breaks spherical symmetry, but a rotating star 
is still axially symmetric. In this case, the TOV equations 
are no longer valid. Hartle and Thorne proposed a per- 
turbative approximation to treat rotating stars (20| . The 
method was further developed and it is valid not only for 
slowly rotating stars, but also for stars rotating up t o the 
Kepler frequency [2l|. In a more recent work [22j the 
Hartle-Thorne approximation was tested against 
a full general relativistic numerical model avail- 
able as part of the numerical relativity library 
LORENE and it was shown that it is reliable for 
most astrophysical applications. In [22] some EoS 
were used to compute rotating star properties, all 
of them obtained at zero temperature. 

In the present work we use the Hartle-Thorne formal- 
ism to calculate the maximum mass, moment of inertia 
and eccentricity of all different classes of possible pulsars 
described above. In previous works [22, [H, [24|, many 
EoS have been investigated, but they were restricted to 
hadronic matter at T = 0. In the present paper we inves- 
tigate all possible classes of pulsars (hadronic, hybrid and 
quark stars) at zero and finite temperature, to account 
also for protoneutron stars. It is important to distinguish 
between the EoS during the short time period when neu- 
trinos are still trapped in the star, and the EoS after the 
neutrinos escape. Next we restrict ourselves to the sec- 
ond case, when pulsars are believed to be already stable 
stars. Notice also that the temperature in the interior 
of the star is not constant @, [25j] > but the entropy per 
baryon is. This is the reason for choosing fixed entropies 
to take the temperature effects into account. The maxi- 
mum entropy per baryon (S) reached in the core of a new 
born star is about 2 (in units of Boltzmann's constant) 
[11. We then use EoS obtained with S = (T = 0), 1 
and 2. 
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Other important motivations for revisiting rotating 
stars is the fact that as they slow down, the decreas- 
ing centrifugal force leads to increasing core pressure and 
density. A softening of the EoS takes place and it can be 
the result of a phase transition to quark matter [27| • If a 
first order phase transition takes place at central densi- 
ties, it was shown that the moments of inertia behave in 
a characteristic way and the braking index diverges [28| . 
Gamma ray bursts are also linked with stellar phase tran- 
sitions and the relation between gamma ray bursts 
and phase transitions in stellar matter under rotation 
should also be extensively considered. The present work 
is the seed for this investigation since we consider only a 
few possibilities for the EoS at T = 0. 

This paper is organized as follows: in Sec. II the for- 
malisms of the Hartle-Thorne approximation are revis- 
ited and the results are presented. In Sec. Ill the results 
are discussed and the main conclusions are drawn. 



II. FORMALISM AND RESULTS 



As the first step we need to know the EoS of the system, 



n = n{p), 



where p is the pressure, e is the energy density, and n 
is the number density of baryons. Once an adequate 
EoS is obtained, it can be used to provide the stellar 
properties. If we opt for a non-rotating configu- 
ration, the Tolman-Oppenheimer-Volkoff (TOV) 
equations [171 ]. which are derived from Einstein's 
equations in the Schwarzschild metric for a static, 
spherical star are used to compute the stellar 
properties, as the radius, gravitational and bary- 
onic masses and central energy densities. If we 
opt for its rotating counterpart, the same EoS is 
used, but the Hartle-Thorne formalism is used 
instead of the TOV approximation. 



where u a is the 4-velocity of the fluid. By using (2) , (3) 
and the Ricci tensor we obtain a system of partial dif- 
ferential equations in which the solution is numerically 
complicated. In the Hartle-Thorne method, rotation is 
treated as a pertubation on the non-rotating configura- 
tion of the star and it gives spherical and quadrupole 
deformations. Within this method the expansion of the 
metric has the form 



e AM) = e A 



1 



to + m 2 P 2 (cos(6')) 



e v [1- 
fl 



2{h + h 2 P 2 (cos(9)))} , 
•2O 2 -/i 2 )P 2 (cos(0))], 



(4) 



where P 2 is the Legendre polynomial of second order; and 
ho, h 2 , too, TO2 and v 2 are all functions of r. 

A non-rotating configuration is computed by integrat- 
ing the Tolman-Oppenheimer-Volkoff equations [l7| of 
hydrostatic equilibrium for the pressure, p(r), and the 
mass interior to a given radius, m(r): 



dp 
dr 
dm 
dr 
dv 
dr 



-(e + p) 



TO + 47IT p 

r(r — 2m) ' 



47rr 2 e, 



(5) 



dr 



where the boundary conditions are to = at r = 0, the 
radius surface R is defined so that v — when r — > 00. 
Here, v is the metric function. In this work we use 
natural units where G = c = 1. 

The quantity U = Q — lu is the angular velocity of the 
fluid relative to the local inertial frame. In the classi- 
cal mechanics the magnitude of the centrifugal force is 
determined by the angular velocity Q. However, in the 
general relativity the magnitude of the centrifugal force 
is determined by 57: 



A. Hartle-Thorne formalism 

The simplest nontrivial form that Einstein's equations 
take is the form for spherically symmetric stars. In this 
case, the line element has the form, 



ds 2 



^)dt 2 + e x ^dr 2 + r 2 (d6 2 +; 



(1) 



However for rotating stars the spherical symmetry is bro- 
ken, but the star maintains axial symmetry. The expres- 
sion for the line element for axially symmetric spacetime 
is given by 



ds 2 = e^dt 2 -e x ^dr 



+ r 2 sin 2 8{d4> - u(r,9)dt 2 ) 2 }. 
The perfect fluid has a stress-energy tensor 
T ah = {e + p)u a u b +pg ab , 



(2) 
(3) 



1 d_ 

r* dr 



dr J 



idj_ 
-—u = 0, 

r dr 



where 



j(r) = e u/2 [l - 2m I r 



1/2 



(G) 



(7) 



The boundary conditions a; = Zu c and dU c /dr = are 
imposed. When r = R we can determine the angular 
velocity, Q, and angular momentum, J, corresponding to 
57 • 



H4f 



n = u(r) 



r=R 



2J_ 
i? 3 ' 



(8) 



To obtain a different value of Vt we need to rescale the 
function ZJ(r): 



(9) 
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The metric function ho obeys the following algebric rela- 
tions 



h 
ho 



mo + J 2 /r 3 



+ 



J 



2m r 3 (r — 2m) 



outside the star, 



-Po + 3 r e 



v u 2 + hoc, inside the star, (10) 



where ho c is a constant determined by imposing ho to be 
continuous at r = R. The mass pertubation factor too 
and the pressure pertubation factor po are calculated by 
integrating: 



G?TO() 

dr 

dpp 
dr 



i de ^ 
dp 

to (1 + 8irr 2 p) 



duj^ 2 



47r(e 



p)r 2 



dr 



(r — 2to) 2 



(r — 2m) 



4 -2 



12 (r - 2m) 



( dm 

\~dr~ 



1 d ( r 3 j 2 u 2 



2m 



3dr 



(11) 



where the boundary conditions are that too(0) = po(0) = 
0. The mass correction at first order is given by 

dm = m + J 2 /R 3 . 

Thus the total mass of a rotating neutron star with cor- 
rections up to the first order is given by 

m(R) +Sm = m{R) + m {R) + J 2 / R 3 , 

where R is the radius of the star surface. 

According to [2fJ, the binding energy of a relativis- 
tic star in a non-rotating configuration is the difference 
between its baryon mass and its total mass-energy 

E B = A - M, 



where A is the total baryon mass 



A = m n / n(r) 
Jo 



'1/2 



4nr dr, 



(12) 



and m n is the nucleon rest mass. To calculate the binding 
energy for the non-rotating configuration we need to solve 
the Tolman-Oppcnhcimcr-Volkoff equation. The first or- 
der correction in the binding energy for rotating config- 
uration is given by 



5E B 
B(r) 



R? 
(e + p)po 

1 



4nr 2 B(r)dr, 

± ( (i - ™) ~ 1/2 - 1 V *i ( i _ ^ ~ 1/2 

dp \ \ r J j dp \ r 



( e ~ ei > ( 1 ~ ~ 



-3/2 



TO.Q 1 2 2—2 

r S 



4irr 2 



1 . 2 4 ( 

— IT 

12 J V dr 



ld f 3-2 

3*™ 



(13) 



where q is the internal energy density, 

6i = e - m n n. 
Thus the baryon mass of rotating stars is given by 

Mb =A + 5E b + 6m. (14) 



The metric functions hi and vi are given by 



where the boundary conditions are /i2(0) = ^2(0) = 0. 
The non-radial mass and pressure perturbation factors, 
m 2 and P2, are determined from the algebric relations 



TO2 = (? — 2m) 
1 



-hi 



1 



dr 



1 



4„-2 



w + -r j 
6 



rZr 



Pi 



u ** 2 — y — 2 



(16) 



To calculate the eccentricity e we use 

1/2 



(17) 



where R p and i? e are the polar and equatorial radii. To 
obtain R p and R e we use the relations 

S(0) = r + £ o (r)+k(r)P 2 (eo S (0)), 
Co = -po{e+p)/(dp/dr), 
£2 = -P2{e+p)/(dp/dr). 



(18) 



where S'(^) is a surface of constant density. 

In table|T]results for slowly rotating stars are presented. 
The M max and R are respectively the maximum mass 
and the radius of an analogous non-rotating configura- 
tion. The M x ax is the maximum mass corrected up to 
first order for a star with angular velocity f2. R is the ra- 
dius of the non-rotating star, R e is the equatorial radius, 
R p is the polar radius, e c is the assumed central energy 
density, / is the calculated moment of inertia and e is the 
eccentricity. The EoS for hadronic and hybrid stars were 
taken from [H and the EoS for quark stars were taken 
from [ll|. As these EoS have already been extensively 
discussed in the literature, we refrain from further ex- 
planations here. The only point worth mentioning refers 
to the parametrizations used. For hadronic and hybrid 
stars we have chosen to work with a parametrization that 
describes the properties of saturating nuclear matter pro- 
posed in [3(| since other common parameter sets for the 
non-linear Walecka model namely, TM1 [31] and NL3 
[32| proved to be inadequate because, due to the inclu- 
sion of hyperons, the nucleon mass becomes negative at 
relatively low densities @, H| ■ The chosen parameters are 
gljml = 11.79 fm 2 , ff 2 /m 2 = 7.148 fm 2 , g 2 p /m 2 p = 4.41 
fm 2 , k/M = 0.005896 and A = -0.0006426, for which 



the binding energy is -16.3 MeV at the saturation density 
po = 0.153 fm^ 1 , the symmetry coefficient is 32.5 MeV, 
the compression modulus is 300 MeV and the effective 
mass is 0.7M. For the meson-hyperon coupling constants 
we have chosen them constrained by the binding of the A 
hyperon in nuclear matter, hypernuclear levels and neu- 
tron star masses (x a = 0.7 and x^ = x p = 0.783) and 
have assumed that the couplin gs to the E and S are equal 
to those of the A hyperon [3fJ, [33| . For the construction 
of the EoS for hybrid stars, the hadronic phase was ob- 
tained with the non-linear Walecka model and the param- 
eters above and the quark phase with the MIT bag model 
with Bag = (180 MeV) 4 . For the quark stars within 
the MIT model, we have used m u = md = 5.5 MeV, 
m s = 150.0 MeV and Bag = (180 MeV) 4 . For the pur- 
pose of the present work, the value of the bag parame- 
ter does not play an important role. Concerning quark 
stars within the NJL model, the set of parameters were 
chosen in order to fit the values in vacuum for the pion 
mass, the pion decay constant, the kaon mass and the 
quark condensates [M A = 631.4 MeV, g s A 2 = 

1.824, g D A 5 = -9.4, m u = m d = 5.6 MeV and m s = 
135.6 MeV which were fitted to the following properties: 
m T = 139 MeV, = 93.0 MeV, m K = 495.7 MeV, 
Ik = 98.9 MeV, (uu) = (dd) = -(246.7 MeV) 3 and 
(ss) = -(266.9 MeV) 3 . 

One can easily see from table U that for all mod- 
els considered, the rotating stars bears a slightly larger 
mass than its non-rotating counterpart. Eccentricities 
are rather uniform for all models (0.43-0.49). Quark stars 
are bound by the strong force and hence they can rotate 
faster than hadronic or hybrid stars that are bound by 
the gravitational force. According to [36j], the value of the 
Kepler frequency can be obtained from the values of the 
mass and the radius of the corresponding non-rotating 
star and its empirical relation reads 



Q = (0.63 - 0.65)(M/ J R : 



3\l/2 



(19) 



In table U the Kepler frequency for each EoS is speci- 
fied. While hadronic and hybrid stars rotate with similar 
frequencies, the values for quark stars are indeed much 
higher. 
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III. RESULTS AND CONCLUSIONS 

Let's now go back to our results in order to compare 
them with what is found in the literature and draw the 
conclusions. 

We first analyzed the results of the rotating pulsars. 
The values shown in table[T]for zero entropy and hadronic 
stars show the same behavior observed in [Hj], i.e., there 
is a small increase in the maximum mass of the rotat- 
ing star as compared with the non-rotating configura- 
tion. Notice, however, that in [2J|, the usual NL3 and 
TM1 parametrizations could be used because only pro- 
tons and neutrons were considered and the central den- 
sity is somewhat lower than in our case. 

In [23[ the models named 0,P,Q and R given in ta- 
bles U can again be compared with our result for the 
hadronic star at T — and they are indeed very sim- 
ilar. The central density and the moment of inertial are 
of the same order and the maximum mass is very simi- 
lar. As far as we know there is no result for rotating stars 
with entropy different from zero (finite temperature) in 
the literature, but we can examine its effect on the stellar 
properties from table |TJ As entropy increases the maxi- 
mum masses decrease for all kinds of pulsars, except for 
the not so precise MIT bag model. Although the radii of 
a non-rotating star degenerally decreases with temper- 
ature (except again for the very simple and unrealistic 
MIT model), the eccentricity of a rotating star remains 
practically unchanged. As expected, the moments of in- 
ertia of quark stars are much lower than of hadronic and 
hybrid stars with consequent larger rotation frequencies. 

The energy released in the conversion mechanism of 
the metastable (MS) star into a stable star (SS) is given 

by 

AE = [(M G {MS) - M G {SS))/M Q ] x 17.88 x 10 53 er 5 . 
where Mq{MS) is the gravitational mass of the 



metastable star and Mq(SS) is the gravitational mass of 
the stable star. For AE to be positive, the gravitational 
mass of the metastable star, at a fixed baryonic mass, has 
to be larger than the gravitational mass of the stable star. 
We have calculated the released energies in the conver- 
sion mechanism for the hadronic to the hybrid, hadronic 
to the quark and hybrid to the quark stars with S = 
in the rotating configuration at fixed baryon mass. The 
released energy is always negative, except in the conver- 
sion mechanism of the hadronic (MS) to the hybrid star 
(SS) at a fixed baryonic mass of 1.56M©, which yields 
AE = 1.14 x 10 51 er<?, allowing for the required energy 
measured in a short gamma ray burst. 

In the present work we have investigated one possible 
variation on the usual static neutral pulsars: the inclu- 
sion of rotation. We have observed that the behaviors 
shown in previous works with much simpler EoS were 
also observed here. The influence of the temperature was 
also investigated in both cases. We are now in a position 
to calculate the energy released from the conversion of 
a metastable star (hadronic or hybrid) to a stable star 
(hybrid or quark) under the influence of slow rotation. 
For the simple calculations done so far with stars in a ro- 
tating configuration at S — 0, only a hadronic star could 
convert into a hybrid star releasing an energy compati- 
ble with a short gamma ray burst. A more detailed and 
complete investigation, with a smaller bag parameter in 
the MIT bag model is under way. 
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TABLE I: Rotating compact stars properties with Kepler angular velocity D. = 0.65(M/7? 3 ) 1 



Type 


Entropy 


M max 


Mi ax 


R 


Re 


R p 


e c 


I 


e 


fi 






(Me) 


(M e ) 


(km) 


(km) 


(km) 


(g/cm 3 ) 


(gem 2 ) 




Hz 


Hadronic 





2.04 


2.07 


11.73 


12.06 


10.86 


1.98 x 10 16 


2.39 x 10 45 


0.43 


1661.37 


Hadronic 


1 


1.96 


2.00 


11.02 


11.33 


10.21 


2.23 x 10 15 


2.11 x 10 45 


0.43 


1791.66 


Hadronic 


2 


1.93 


1.96 


10.91 


11.22 


10.11 


2.24 x 10 15 


2.03 x 10 45 


0.43 


1801.74 


Hybrid 





1.63 


1.65 


12.33 


12.78 


11.19 


1.57 x 10 15 


1.90 x 10 45 


0.48 


1381.09 


Hybrid 


1 


1.50 


1.53 


11.32 


11.73 


10.28 


1.75 x 10 15 


1.56 x 10 45 


0.48 


1503.07 


Hybrid 


2 


1.50 


1.52 


11.76 


12.20 


10.65 


1.58 x 10 15 


1.65 x 10 45 


0.49 


1418.80 


v^uarKomc^ivii i ) 


u 


1 00 
l.ZZ 


1 0^ 
1 .ZD 


0. / 




0*3 


0.14 X 1U 


U.Oo X 1U 


U.40 


OQQO QQ 

Zyoo.oy 


Quarkonic(MIT) 


i 


1.22 


1.25 


6.76 


6.97 


6.23 


5.17 x 10 15 


0.53 x 10 45 


0.45 


2941.47 


Quarkonic(MIT) 


2 


1.23 


1.26 


6.79 


6.99 


6.25 


5.08 x 10 15 


0.54 x 10 45 


0.45 


2929.40 


Quarkonic(NJL) 





1.20 


1.23 


7.87 


8.12 


7.24 


3.45 x 10 15 


0.69 x 10 45 


0.45 


2316.71 


Quarkonic(NJL) 


1 


1.17 


1.18 


7.71 


7.81 


7.45 


3.68 x 10 15 


0.64 x 10 45 


0.45 


2363.51 


Quarkonic(NJL) 


2 


1.10 


1.13 


7.18 


7.39 


6.62 


4.58 x 10 15 


0.51 x 10 45 


0.44 


2552.36 



